Geometric Stability Spectra of Dipolar Bose Gases in Tunable Optical Lattices 
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We examine the stabiHty of quasi-two-dimensional dipolar Bose-Einstein condensates in the pres- 
ence of weak optical lattices of various geometries. We find that when the condensate possesses a 
roton-maxon quasiparticle dispersion, the conditions for stability exhibit a strong dependence both 
on the lattice geometry and the polarization tilt. This results in rich structures in the system's 
stability diagram akin to spectroscopic signatures. We show how these structures originate from the 
mode matching of rotons to the perturbing lattice. In the case of a one-dimensional lattice, some 
of the features emerge only when the polarization axis is tilted into the plane of the condensate. 
Our results suggest that the stability diagram may be used as a novel means to spectroscopically 
measure rotons in dipolar condensates. 

PACS numbers: 67.85.Hj, 03.75.Hh , 05.30.-d 
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I. INTRODUCTION 

The field of ultracold atoms has become a useful prov- 
ing ground for theories of condensed matter systems. 
Progress with Bose-Einstein condensation in particular, 
both experimental and theoretical, has increased our un- 
derstanding of superfluidity and long-range order and 
their dependence on interactions [l|. Recent years have 
seen increasing attention given to condensed systems 
whose constituent bosons interact via the dipole-dipole 
force, which is long-range and anisotropic. Such dipolar 
Bose-Einstein condensates (dBECs) have been achieved 
with atomic ^^Cr 0, ^^^Er [sj, and ^^"^Dy [4], and re- 
searchers are making active progress toward the conden- 
sations of more strongly interacting polar molecules [5|- 
0] and Rydberg atoms [8|. The atomic dBEC experi- 
ments have observed various consequences of the dipole- 
dipole interaction, includinggeometry-dependent stabil- 
ity |3, Q, c^-wave collapse [sl, UqI: ^^^ deconfinement- 
induced collapse [11]. 

Dipolar condensates are predicted to exhibit a roton- 
maxon structure in their dispersion relations p^2l-[l^, 
qualitatively similar to that observed in superfiuid ^He 
half a century ago [15l-[l8j. In the case of dBECs, the 
depth and location of the roton minimum can be tuned 
by variations of the interaction strength, density, and 
trapping parameters. A system becomes dynamically un- 
stable when the roton energy is tuned to zero, although 
quantum and thermal fluctuations may destroy the con- 
densate for small, but finite, roton energies [i9|,[20|. Var- 
ious signatures of rotonization have been predicted re- 
cently, such as reduced critical superfiuid velocity [12], 
abrupt transitions in Faraday patterns |2l| , ro ton-length- 
scale oscillations of two-body correlations [22] , roton con- 
finement [23], and short- wavelength immiscibility phases 
|2^ . To date, none of these effects have been observed 
in dBECs, although a variant of Bragg spectroscopy 
recently measured the roton of a nondipolar gas with 
cavity- mediated long-range interactions [25[ . Bragg spec- 
troscopy, as calculated in Ref. J26|, should therefore be 



feasible in dipolar condensates. 

The recently proposed method of stability spec- 
troscopy [27] is an appealing alternative for roton mea- 
surement that produces multiple signatures, while ex- 
ploiting the inherent instability and sensitivity of ro- 
tonized systems. The general idea is to probe the ro- 
tonized dBEC with a weak lattice of periodicity Al and 
depth s, and then straightforwardly map the stability as 
a function of these lattice parameters. For quasi-one- 
dimensional (qlD) dBECs, it was demonstrated that the 
critical lattice depth 5c, above which no stable conden- 
sate exists, depends strongly on the periodicity and ex- 
hibits local minima whenever Al equals either half the 
roton wavelength or a low roton subharmonic. The sta- 
bility boundary 5c (Al) hence constitutes a spectroscopic 
measurement of the roton wavelength, and it is easily 
measured because the lattice depth and spacing can be 
adjusted by varying respectively the intensity and angle 
of crossed off-resonant laser beams [28]. We emphasize 
that the necessary stability measurement is effectively 
binary: the condensate is either stable or unstable. 

The goal of the present article is to generalize and ex- 
tend the method of stability spectroscopy to quasi-two- 
dimensional (q2D) dBECs, elucidating the crucial roles 
played by lattice geometry and polarization tilt. We will 
show that, in a ID lattice, the central feature of the qlD 
case (at Al = Arot) persists in q2D, but the other fea- 
tures emerge only when the polarization axis is tilted 
into the trapping plane. Additionally, we consider the 
stability spectra that result from 2D triangular-lattice 
perturbations, where multiple roton signatures are most 
pronounced when the polarization is orthogonal to the 
trapping plane. Although our results are computed nu- 
merically, we explain the origins of the relevant stability 
structures by applying perturbation theory to the Gross- 
Pitaevskii and Bogoliubov de Gennes equations. Section 
nil describes our mean-field formalism. Sec. IIIII outlines 
the perturbation theory employed in our analysis. Sec. 
IIVI presents and explains the stability spectra for several 
interesting q2D cases, and Sec. |Vl concludes our discus- 



sion. 



II. MEAN-FIELD FORMALISM 



with a and 77 being respectively the polar and azimuthal 
angles defining the polarization direction d ^]. Under 
the single-mode approximation, the ground state satisfies 
the q2D Gross-Pitaevskii equation 



We consider a dilute gas of N interacting bosons that 
is tightly confined in the z direction by a harmonic trap 
of frequency ujt and moves freely in the xy plane. The 
system is then subjected to a weak lattice perturbation in 
the xy plane represented by U{p). We scale lengths and 
energies by the natural units of the problem, which are 
It = y/h/mujt and HuJt respectively. In the mean-field 
description of the system, the perturbed ground state 
^o(r) minimizes the grand canonical Hamiltonian 



Hm = /(iV^*(r) 



1 1 



-^N Jd'r'l^ir'fVsDir-r') 



(1) 



^(r) 



where /j^sd is the chemical potential and V3D(r) is the 
interaction pseudopotential for polarized dipoles [29j 



^3D(r) = 4.7rasS (r) + 3add- 



{d-ry 



(2) 



The interaction potential depends on the scattering 
length as and the dipole length add = md^ /?>h?^ as well 
as the polarization direction d. The function ^(r) is as- 
sumed to have unit norm. 

We make the simplifying assumption that both U{p) 
and the interaction energy of the particles are small com- 
pared to the tight trapping energy. This allows for the 
single- mode approximation [30], in which the order pa- 
rameter is confined to zero-point oscillations along the 
z direction. The z dependence of the order parameter 
then factorizes out as ^(r) = Ti~^e~^ ^'^ip{p). With this 
ansatz, Eq. ([!]) reduces to 



H[^] 



d^prip) 



\K + u{p) 
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+\n I d^p'\^l^{p')\^V{p-p') 
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where the q2D interaction V is given in momentum space 
as 



y(k) 



47r 



as + add cos aF±_ 



V2) 



2 / k \ 

-sin aF\i I — I 



V2J 



(4) 



The functions F±{q) and -F'ii(q) are defined by 
-F±(q) = 2 - 3\/7rge«'erfc(g) 
F\\{^) = -1 + 3v7r^ -e"^ erfc{g) 



(5) 



^^Mp) = - 7:KMp) + u{p)Mp) 



+ N d^p'\Mp')\ vip-p')Mp) 



(6) 



with the effective q2D chemical potential /i = /isd — \' 

The dynamical excitations above the ground state ijjq 
are determined by solving the Bogoliubov de Gennes 
equations. These can be written compactly as 



H. 
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E, 



J 
(7) 
where Hq = ~^V^ + U{p) is the noninteracting single- 
particle Hamiltonian, C[x]{p) = N j d^ p' |V^o(pOI^ V^P- 
p')x{p) describes direct interactions with the condensate, 
and X[x]{p) =Njd^p'x{p')Mp')V{p - p')Mp) is an 
integral operator describing exchange interactions. The 
functions Uj and Vj are subject to the usual normalization 

condition J d'^p (\uj\ — \vj\ ) = 1- The ground state i/jq 

satisfying Eq. ([6]) is dynamically unstable if one or more 
of the excitation energies Ej is imaginary- valued, causing 
local collapse on a length scale set by the unstable mode 

In the absence of the perturbation U{p)^ the solu- 
tions to Eq.® and Eq. ([7|) are well known. Transla- 
tional invariance guarantees that momentum is a good 
quantum number. For simplicity, we will assume peri- 
odic boundary conditions over a rectangular domain of 
area A, which discretizes the momenta of the system. 
It is also convenient to introduce the integrated (over 
z) density n2D = N/A. The unperturbed ground state 
and chemical potential are then '^q \p) = 1/a/A and 



p 



(0) 



n2DV{k) 



k=0 



respectively. The excitations are 

parameterized by their momentum quantum number k, 
and are given by [33| 
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where (p\^{p) = e*^"^/V^ and E^ 



(0) 



(8) 



c.(k) 



k'^/2 lk'^/2-\- 2n2DV(k)] is the Bogoliubov spectrum. 

For certain densities and interaction parameters, the 
dispersion relation uj{k.) may contain a local minimum 
at nonzero momentum. The corresponding mode is re- 
ferred to as a roton mode. When the dipoles are polar- 
ized orthogonal to the plane of motion, the interaction 
and dispersion depend only on the momentum magni- 
tude /c, thereby causing the set of roton modes to form a 
ring of radius k^ot in /c-space. Tilting the dipoles into the 



plane, however, makes both the interaction and disperion 
anisotropic. As an example, consider the case in which 
the polarization is tilted somewhat towards the y axis. 
Modes that propagate along y create density antinodes 
along lines of constant x, thereby accumulating dipoles in 
the higher-energy side-to-side configuration. Conversely, 
modes that propagate along x tend to accumulate dipoles 
along lines of constant y, where the dipoles are some- 
what head-to-tail because of the tilt. In a general sense, 
modes that propagate along the tilt projection tend to 
have higher energy than those that propagate in a di- 
rection orthogonal to the tilt projection, and there is a 
continuous transition between the two as the direction 
of a mode is varied. This anisotropy of dispersion has 
been shown to lead to anisotropic superfiuidity [31] and 
coherence [34], as well as striped immiscibility states in 
binary dipolar condensates [2J]. We will see in Section 
IIVAI that this anisotropy leads to emergent features in 
the stability spectrum of a q2D dipolar condensate in a 
ID lattice. 



III. PERTURBATION THEORY 

Reference [27] demonstrates that the main features 
of a qlD stability spectrum can be understood in the 
context of simple perturbation theory. We have found 
that this is also true for q2D spectra. In this section, 
we briefly develop the essentials of perturbation theory 
which are necessary to understand the results of Sec. 
llVi We will establish first-order perturbative results for 
the order parameter and mean-field potential, and then 
write down the first-order perturbation equation for Bo- 
goliubov modes. Further details regarding rigorous per- 
turbation theories applied to the Gross-Pitaevskii and 
Bogohubov de Gennes equations can be found in Refs. 

To study the response of the condensate to the lat- 
tice perturbation /7(p), we expand both the order pa- 
rameter and chemical potential in perturbation series as 
^ = ^(0) _|_ ^(1) ^ _ , and II = /i^^) + /i^^) + . . . , respec- 
tively. To first order, the Gross-Pitaevskii equation (|6j) 
reduces to the linear perturbation equation 



,(1)^(0) =_lv>(^)(p) + [/(p)V^(«)(p) 



2n 



ID 



j cPp'i^^'\p')V{ 



P-P') 



(9) 



under the requirement that the order parameter remains 
normalized to first order, implying that ('^^^^l'^^^^) = 0. 
Multiplying both sides of Eq. ([9]) by ip^^^ and integrating 
yields the correction to the chemical potential 



/,(!) = /^(0)|C/|^(0) 



(10) 



We compute the correction to the order parameter by 
expanding ip^^^ in the complete basis of complex expo- 
nential functions {(^\<i} introduced below Eq. (|8|). We 



then find that the expansion amplitudes are given by 



(^kl'*/'^ 



£(k) 



(11) 



where £(k) = /c^/2 + 2niDV(k) is the Hartree-Fock en- 
ergy [40]. Given that cj(k) = ^^e{k)k'^ /2^ the presence 
of a roton mode in the Bogoliubov dispersion implies a 
local minimum in the Hartree-Fock spectrum at a com- 
parable value of k. For soft-roton systems, the locations 
of these minima very nearly coincide. It is easy to show 
that £:(k) -^ 2/i^^^ in the zero-momentum limit, indicat- 
ing that Eq. (pT]) is well defined even for small k when- 
ever the unperturbed condensate is stable. 

For the purposes of this paper, it is useful to define 
the combined potential experienced by an atom at posi- 
tion p, which is the sum of the external and mean-field 
potentials. We evaluate this potential to first order from 
([TT]) . which results in 



Ue{p)^U{p)-n + Njd^p'mp')fv{p-p') 



kT^O 



kp_ 



2£(k) 



^k\U\^^°A + O {U 



(12) 



after some algebra. A typical lattice potential projects 
onto only a few basis modes (/^k, which greatly simpli- 
fies the formulae for both '^^^^ and Uc in practice. In 
the absence of interactions, the combined potential is of 
course equal to the external potential itself. However, 
dipolar interactions may cause the system to rotonize, in 
which case the mean-field potential may amplify the per- 
turbation by a factor of order /c^/2£(k), the magnitude 
of which depends strongly on the modes that compose 
the lattice. 

To study the response of quasiparticle energies to the 
perturbation, we now apply a simple perturbation theory 
to the Bogoliubov de Gennes equations. We expand Ej, 
Uj^ and Vj in perturbation series, and then substitute into 
Eq. ([7]) to find the perturbation equation: 









/7i'^+X(i) X(i) 

x(i) [/P+x(i) 





^(1) f 4°^ 



where X^^^[x]{p) = N J £p'x{p')^I^^^^^V{p - p') 
describes exchange interactions with the unper- 
turbed condensate, Uc (p) is the first-order combined 
potential given by Eq. ([Ej), and X^^^[x](p) = 
Njd^p'x{p')i^^°^V{p - p'){i;<-'Hp')+i;<-'Hp)) de- 
scribes exchange interactions with the condensate 
perturbation. Our present interest lies only in the 
energy shift E^^ , which we can isolate by acting on both 
sides of Eq. (p!3|) by the operator J d^p (u\^f ^v^./ J for 
any k' satisfying the degeneracy condition uj{k.^) = cj(k). 




FIG. 1: (Color online) Stability spectra for q2D dBECs in a weak lattice with varied polarization tilts and scattering lengths, 
chosen in such a way that uj(kx, 0) is similar in each case. The shaded regions correspond to stable, condensed ground states. 
The experimental parameters are for ^^^Dy trapped at cjt = Stt x 10 kHz with a density nsn ~ 10^^ cm~^ and dipole 
moment d = 10/is. Cases I-IV correspond respectively to polarization tilt a (towards y) and scattering length as values 
oi {a,as}i = {0°,— 82ao}, {a,as}ii = {7°, — 79.7ao}, {a,as}ui = {10°, — 77.35ao}, and {a,as}iv = {15°, — 71.6ao}. The 
right-hand panels depict the dispersion relations for these cases, with the blue (solid) lines corresponding to (jj{kx,0) and the 
red (dashed) lines corresponding to uj{0,ky). For each case, the roton wavelength (for modes parallel to x) is approximately 
Arot ~ 575 nm. The vertical dashed lines on the stability spectrum denote Arot/2, Arot, and 2Arot, which approximately identify 
the spectral features. 



First-order energy shifts are then determined by the 
Hermitian matrix on the right-hand side of Eq. (p!3|) , 
which we denote compactly as 



A: 



U, 



(1) 



• x(i) 



xd) 




(14) 



Our analysis will focus primarily on the softening of roton 
modes, which generally are closely connected to system 
instability. 



IV. STABILITY SPECTRA 

A. One-Dimensional Lattice 

We now consider the stability spectrum of a rotonized 
q2D dipolar condensate that is perturbed by a one- 
dimensional lattice of tunable depth and lattice spacing. 
Without loss of generality, we focus on the special case 
in which the lattice is directed along the x-axis: 



U{p) = s cos {ki,x) . 



(15) 



Such a potential may be generated by a pair of crossed 
off-resonant beams of wavelength ALas and angle 6>, with 
the lattice spacing determined by the relation Al = 
27r//cL = ALas/2sin(6>/2). The depth parameter is 
proportional to the single-beam intensity Iq via s = 



—Re{a{uj)}Io/4:eochuJt^ where Re{a{uj)} is the atomic 
polarizability. The depth and spacing of the lattice may 
be tuned by respectively varying the intensity /q and an- 
gle [23. For each spacing Al, the dBEC destabilizes at 
depths above a certain critical value 5crit(AL). 

Our numerical method maps the stability of a system 
by exploiting two intrinsic symmetries of the ID lattice 
(p!5|) . First of all, the system remains invariant with re- 
spect to translations in y. The order parameter of a 
stable perturbed ground state should thus be indepen- 
dent of this variable, although a collapsed state will gen- 
erally exhibit a spontaneous breaking of this symmetry 
due to local collapse [32[- With this in mind, we nu- 
merically solve the Gross-Pitaevskii equation (|6]) using 
the conjugate gradients method [41] under the assump- 
tion of translational invar iance along y. We then use 
the computed ground state to numerically solve the Bo- 
goliubov de Gennes equations ([7j), which themselves are 
block-diagonal due to both parity in x and translational 
invar iance in y. The signature of a collapsed state is the 
existence of imaginary- valued Bogoliubov excitation en- 
ergies for an assumed ^/-independent order parameter, or 
alternatively local collapse. Our numerical grid consists 
of 2^ points along each of x^y e [— 24Arot, 24Arot]- 

Figure [1] shows the numerically-computed stability 
spectra for several different rotonized dBECs. The dif- 
ferent stability boundaries correspond to different polar- 
ization tilts towards the y-ax.is (77 = 90°, in Eq. ([5|)), 
whose dispersions uj{kx^O) and uj{0^ky) are depicted in 



the subplots. For each case, we have chosen the tilt ol 
and scattering length a^ in such a way that the disper- 
sions along X, most particularly the depth of the roton 
mode, are comparable in each case as shown by the solid 
blue lines. This allows us to isolate the consequences of 
polarization tilt from those related to roton depth. As 
described earlier, polarization tilt causes the modes di- 
rected along the tilt projection (^ in this case) to have 
higher energy than modes that are orthogonal to the tilt 
projection (along x). 

In Fig. [H we observe that the central stability feature 
at Al ~ Arot appears to exist independent of polarization 
tilt. It originates in the strong amplification of the mean- 
field potential, as occurs in qlD systems [27]. This may 
be seen by evaluating the combined potential ([12]) given 
the ID lattice (fTSj) : 



Uc{p) 



yleiki^x) 



kl 



cos(A:lx) + (9(s^). 



(16) 



In each case considered, the Hartree-Fock energy £:(k) 
has a shallow local minimum at ki,x ~ k^ot^ due to ro- 
tonization. Equation (p!6|) thus indicates that the mean- 
field strongly amplifies the perturbation by a factor of 
k'^^^/2£{krot^) ^ 1 for Al ~ Arot- This mean-field en- 
hancement of the lattice probe is responsible for the cen- 
tral stability dip in each case. A similar effect was re- 
cently observed in simulations of nonrotonized dBEC su- 
perfiow in a weak lattice [42]. In that context, mean- 
field enhancement originates from the polarization tilt 
(towards y) effectively lowering the trapped-gas analogue 
of eiki^x) for fixed scattering length and lattice spacing. 

The appearance of stability dips at Al = Arot/2 and 
Al ~ 2Arot depends strongly on the polarization tilt an- 
gle a. This can be understood in terms of the mode- 
matching of Bogoliubov roton modes by the Hermitian 

matrix A defined in Eq. ^^. Both /(/:'k'|^i^^l^k) 

and {(p^^f\X^^'^\(p^^) are proportional to 4;,fe^^|/c^-fe,|,/cL, 
which defines the mode-matching conditions of the per- 
turbation. Degenerate modes satisfying these conditions 
experience a first-order energy shift. Given the reflection 
symmetry of the perturbation (p!5]) . a mode-matched pair 
will separate into even and odd solutions in x. For s > 0, 
the odd solution will lower in energy because it accumu- 
lates atoms in the minima of the combined potential (p!6|) . 
as illustrated in Ref. [27]; conversely, the even solution 
will increase in energy. This mixing and splitting of de- 
generate, matched modes is similar to that of "staggered 
modes" which appear at the Brillouin zone edge in single- 
particle band theory [43] , although we are considering the 
dispersion of quasiparticles. 

Because roton modes are local minima of the system 
dispersion, we expect hastened destabilization for all lat- 
tice periodicities in which low-energy roton modes soften 
to first order. Figure [2] depicts the /c-space ring of degen- 
erate rotons for the a = case. For /cl G (0, 2/crot), there 




FIG. 2: (Color online) Depiction of roton modes for the case 
of zero polarization tilt (a = 0). Degenerate modes satisfying 
ky = k'y and \kx — k'x\ = ki, (indicated by pairs of green dots) 
experience first-order energy shifts. 



are always two pairs of rotons. 



±A:L/2,Vfc?ot-(fcL/2r 



and 



±fcL/2, ■ 



^rot 



(fcL/2)' 



which mode match. For /cl = 2/crot, these two pairs 
coalesce into a single pair at k = (±/crot,0). The fact 
that there is always a low-energy roton pair to soften for 
Al ^ Arot/2 implies that the BEG is fairly unstable to 
collapse at small s for all applied lattice wavelengths in 
this range. The stability boundary I in Fig. [T]is therefore 
low and relatively featureless except for the central dip 
already accounted for. By contrast, tilting the polariza- 
tion towards y raises the relative energies of rotons with 
a nonzero ky component. With increasing tilt a (and 
fixed cj(/crot, 0)), we then expect an emergent stability dip 
when the lowest-eneigy rotons (along ±f ) mode-match. 
This occurs when k\, = 2krot (Al = Arot/2), explaining 
the other prominent feature in Fig. [1] case IV. This is 
in direct analogy with the corresponding stability dip in 
the qlD scenario [23, in which case excitations with a 
y component have higher relative energies due to tight 
trapping, rather than polarization tilt. We note that we 
have observed the emergence of the Arot/2 stability fea- 
ture in fully-3D simulations of fiat pancake dBEGs, where 
Arot is less well defined. Thus, the emergent stability dip 
is observable in the experimentally-realistic scenario of a 
rotonized trapped gas. 

The smaller feature at Al ^ 2Arot in Fig. [1] derives 
from the interplay of second-order effects, such as (mod- 
est) mean-field amplification and the mode-matching of 
degenerate rotons and phonons along x. The shift to 
a slightly shorter wavelength is a consequence of the 
Hartree-Fock minimum occuring at a somewhat larger 
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FIG. 3: (Color online) Schematic for the alignment of lasers 
(in the xy plane) used to produce the triangular lattice that 
we consider. 



momentum than /Crot- Moreover, the first mode mixture 
to destabilize in this case may not even include the pre- 
cise roton minimum, since the (degenerate) second-order 
energy shifts responsible for destabilization depend on 
the dispersion as a whole. This complicates the precise 
analytical determination of this feature's location beyond 
the approximation Al ^ 2Arot; however, this approxima- 
tion improves as the unperturbed roton mode softens in 
energy. 

The stability structures discussed in this section tend 
to disappear when the polarization is tilted towards the 
lattice vector (along x). Such a tilt increases the Hartree- 
Fock energy e (k\,x\ thereby reducing the effect of mean- 
field amplification in Eq. (p!6|) . Moreover, the lowest- 
energy rotons (now along ±^) only mode-match to each 
other when k\, is vanishingly small. As a result, the sta- 
bility boundary of such cases tends to be relatively struc- 
tureless and is not depicted here. 



B. Triangular Lattice 

For the case of a ID lattice imposed on a rotonized 
q2D dBEC, a polarization tilt is required for the appear- 
ance of multiple stability dips. We now demonstrate that 
a triangular lattice perturbation causes multiple stabil- 
ity structures to emerge even in the a = case. One 
may produce a triangular lattice using three crossed, off- 
resonant beams as depicted in Fig. [3l Without loss of 
generality, we assume a stationary beam directed along 
X, with symmetrically-placed beams making an incidence 
angle 7 with the x-axis as shown [44]. Assuming equal 
polarization and intensity among the beams, the corre- 
sponding potential is [45] 

[/(p) = s [2 cos (/cLas(l + cos7)x) COS (A:Las(sin 7)?/) 
+ cos(2A:Las(sin7)?/)] 

(17) 
where /cLas = 27r/ALas is the laser wave number and s is 
the AC Stark shift factor defined below Eq. ^ [46]. 
As the notion of lattice spacing is ambiguous for this 
potential, we instead vary the crossing angle 7 G (0, 180)° 
and identify the stability boundary Scrit(7)- 

Due to the lack of continuous translational invariance 
of the triangular lattice ([T7|), the block-diagonals of the 
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FIG. 4: (Color online) Stability Spectrum for a (zero- 
polarization-tilt) q2D dBEC subjected to a triangular lattice 
perturbation. The shaded region corresponds to a stable, con- 
densed ground state. The experimental parameters are iden- 
tical to those of case I in Fig. [U The assumed laser wave- 
length is ALas = 595 nm. The vertical dotted lines denote 
the predicted locations of stability dips, based on first-order 
perturbation theory applied to the Gross Pitaevskii equation. 
We have employed piecewise cubic interpolation to smooth 
the stability boundary. 



Bogoluibov de Gennes eigenvalue problem ([7]) are large 
compared to those of the ID lattice scenario in Sec. IIVAI 
This makes repeated diagonalization more computation- 
ally expensive. We instead assess stability directly from 
the Gross-Pitaevskii energy functional (J3j). Implement- 
ing a conjugate gradients algorithm [4l|, we compute the 
energy minimum (converged to a relative error of 10~^^, 
if a minimum exists) and use local collapse as the sig- 
nature of instability. Our simulations use a grid of 2^ 
points along x G [— IGArot, 16Arot] and 2^ points along 
y G [— SArot, SArot]- We choose a fixed, horizontal laser 
mode /cLas from our /c-space grid, and compute stabil- 
ity only for angles 7 associated with /c-grid modes sat- 
isfying ||k|-A:Las| < dk:j,l2, where dk^^ = 7r/16Arot is 
the numerical grid spacing of kx modes. This enforces 
in an approximate way our assumption that the three 
lasers forming the lattice have equivalent wavelengths, 
while guaranteeing that the lattice satisfies the assumed 
periodic boundary conditions of our grid. We further 
employ shape-preserving piecewise cubic interpolation to 
the smoothly-varying stability boundary for angles that 
cannot be produced on our finite, periodic spatial grid. 
We find that these numerical approximations are suffi- 
cient to resolve the important spectral features of the 
stability diagram. 

Figure H] shows the triangle-lattice stability spectrum 
of the same a = dBEC as was used in Fig. [H We 
observe three main features, each of which originates en- 
tirely from mean-field amplification. We compute the 
first-order combined potential by substituting Eq. (pTj) 




FIG. 5: (Color online) Predicted locations of stability struc- 
tures as a function of laser wavelength ALas- The symmetric 
blue lines are solutions to sin 7 = ALas/2Arot, and the red 
line is the solution to cos ^ = ALas/2Arot. Angles correspond- 
ing to the laser wavelength used for Fig. 3] are indicated by 
green dots, and they define the locations of the vertical lines 
overlain on the triangle-lattice stability spectrum in Fig. 2] 



into Eq. ([I2|): 

Uc{p) = s [2A/eL,,,7 COS (/cLas(l + cos7)x) COS (A:Las(sm 7)^) 

+^feLas,7Cos(2A:Las(sm7)?/)] + O {s^) 

(18) 
where 

Ai 

^fcLas,7 

(19) 



B, 



fcLas,7 



(2/CLasCOs(^))' 
2s(2fcLasCOs(^)) 

(2feLasSin7)^ 

2£(2/cLasSin7) 



Given that s{k) has a shallow local minimum at /c ~ /crot, 
a component of the perturbation (pTj) is strongly ampli- 
fied when either sin 7 ^ A:rot/2A:Las or cos(^) ^ A:rot/2A:Las 
is satisfied. Figure [5] plots the three solutions to these 
equations on the domain 7 G (0, 180)° as functions of 

ALas/Arot = /^rot/^^Las, which Cxist Ouly for ALas < 2Arot 

(or, equivalent ly, /cLas > A:rot/2). For the particular laser 
wavelength relevant to Fig. H] (ALas/Arot = 1.03), three 
solutions are identified (green dots) which determine the 
values of 7 where stability minima should occur. These 
predicted minima are depicted in Fig. 2] with vertical 
lines. As we can see, the observed locations of the stabil- 
ity dips agree quite well with the prediction. The stabil- 
ity boundary near 7 ^ 120° is somewhat lower because 
it originates from the mean- field enhancement of the first 
term in Eq. ([T8|) , which has an additional factor of 2. 

It is interesting that these locations can be tuned by 
varying the laser wavelength, with the constraint that the 
wavelength remain detuned from internal atomic transi- 
tions. We note that the two structures which are sym- 
metric with respect to the angle 90° (located at 7s < 90° 



and 180° — 7s > 90°) would appear even in the absence 
of the horizontal beam, in which case the perturbation 
reduces to a ID lattice with a periodicity that is sym- 
metric over 7 = 90°. The asymmetric structure (whose 
location we denote by 7^) would appear in the absence 
of either of the oblique laser beams for similar reasons. 
The simulataneous appearance of all three structures re- 
quires the 2D character of the lattice. Once the loca- 
tions of these stability dips are measured, one can infer 
the roton wavelength via either Arot ^ ALas/ (2 sin 7s) or 

Arot ^ ALas/ (2 COS ^) . 

Similar to the ID-lattice scenario in Sec. IIV A( mode- 
matching plays a minimal role when a zero-tilt rotonized 
dBEC is perturbed by a triangular lattice. This may 
be understood by considering matrix elements of the A 
matrix (defined below Eq. (p!3|) ) between modes k' and k. 

In this case, both ((p\i'\Uc |^k) and ((/?kH-^*^"^^l^k) are 
proportional to a sum of the Kronecker delta functions 



and 



-^\k'^ -kx I ,/cLas (1+cos j)0\k' -ky I ,/cLas sin 7 



^fc^,/c^^|/cJ^-/c^|,2/cLas sin77 



(20) 



(21) 



where mode-matching can occur essentially via either 
term in Eq. (p!8|) . This is in direct analogy with Fig. 
[2] except that matched modes connect either vertically 
by the vector ±2/cLasSin7^ or diagonally by the vec- 
tor ±/cLas(l + cos7)f ± /cLas sin7^, the norm of which is 
2A:Las cos ^. Rotons cannot mode-match if both of these 
vectors are too large in magnitude to form chords in a 
circle of diameter 2/crot- Such is the case for angles in 
the vicinity of 7 = 60° when /cLas > 2/crot/\/3 (equiva- 
lently, ALas < A/3Arot/2), as follows from geometric prin- 
ciples. This raises the stability boundary somewhat for 
such angles, since the destabilization of rotons is then 
a second-order effect; nevertheless, no new stability fea- 
tures emerge. 

The question naturally arises as to how polarization tilt 
might affect the stability spectrum of a rotonized gas in 
a triangular lattice. We have found that certain features 
in Fig. Slwash out with increasing tilt, depending both 
on AL/Arot and the direction of polarization tilt. For the 
case of tilted dipoles, the amplification factors in Eq. (p!9|) 
are written more generally as 



^A;Las,7 ~ 



2e{ 



(2/CLasCOs(^)) 

(±A:Las(l + COS7)x ± /CLas siu 

_ (2feLasSin7)^ 
^^^''^~ 2£(±2^LasSin7^) 



ly) 



(22) 



Since the polarization tilt raises the relative Hartree-Fock 
energies of roton modes along the tilt projection, this re- 
sults in A/cLas,7 and/or B^^^^^^ being less-strongly-peaked 
functions of 7. For example, tilting the polarization to- 
wards y raises the relative value of e{±kroty), which di- 
minishes the enhancement due to ^/cLas 7- This would 



cause a washing out of the symmetricany-located fea- 
tures in Fig. m Similarly, a polarization tilt towards x 
would diminish the enhancement due to A/cLas,7 ^^^ wash 
out the central stability feature. 

We note that the stability structures originating from 
mean-field enhancement, both in the ID and 2D lattices, 
can be observed in systems whose dispersions are not 
quite rotonized (ie, possess no roton minimum in co'(k)). 
This is because the Hartree-Fock energy £(k), whose in- 
verse appears in the combined potentials ([16]) and (p!8|) . 
may have a local minimum even when the dispersion 
co'(k) does not exhibit a roton. In such cases, the in- 
verse Hartree-Fock energy l/^(k) is less strongly peaked, 
resulting in respective stability dips that are generally ob- 
servable, albeit less pronounced. The locations of these 
dips are set by the local minimum of 6:(k), and this is 
strictly true even for rotonized systems, in which case krot 
only approximates the location of the Hartree-Fock min- 
imum. These particular structures thus do not measure 
the onset of rotonization in the dispersion, but rather the 
closely-related phenomenon of a Hartree-Fock local min- 
imum at nonzero momentum. Such can occur at lower 
densities than are generally necessary for rotonization, 
making our results directly relevant to current dBEC ex- 
periments that do not appear to have fully rotonized dis- 
persions. 



tonized dipolar condensate is perturbed by a weak optical 
lattice, the resulting stability plot contains spectroscopic 
information from which the roton wavelength may be in- 
ferred. In this paper, we have examined the stability 
spectra of rotonized q2D dipolar BECs in the presence of 
tunable ID and 2D lattices. The tilt of the dipole polar- 
ization axis plays an important role in determining which 
structures are present in the ID-lattice stability diagram. 
The central stability feature at periodicity Al ~ Arot ex- 
ists independently of polarization tilt, whereas the fea- 
tures at Arot/2 and 2 Arot emerge only when the rotonized 
gas has the polarization axis tilted into the plane, re- 
maining orthogonal to the lattice vector. The emergent 
structures originate in the mode-matching of directional 
rotons, the relative energies of which depend on polar- 
ization tilt. In the case of a triangular lattice perturb- 
ing a zero-tilt dBEC, we find structure in the stability 
spectrum resulting from the 2D character of the pertur- 
bation. The locations of stability dips can be tuned by 
varying the off-resonant laser wavelength used to pro- 
duce the lattice. Recent advances in the condensation of 
magnetically dipolar atoms suggest that measurements 
of stability spectra should be feasible, and that rotons in 
dBECs may be consequently observed. 
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